Preprint typeset in JHEP style - HYPER VERSION 



HU-EP-07/29 



Anomalous dimensions of finite size field strength 
operators in J\f = 4 SYM 



Matteo Beccaria 

Dipartimento di Fisica, Universita' del Salento, Via Arnesano, 73100 Lecce 
INFN, Sezione di Lecce 



E-mail: matteo .beccaria@le . infn . it 



Valentina Forini 

Humboldt-Universitdt zu Berlin, Institutfiir Physik, Newtonstrafie 15, D-12489 Berlin 
and 

Dipartimento di Fisica, Universita' di Perugia, Via A. Pascoli, 1-06123 Perugia 
INFN, Sezione di Perugia 
E-mail: 



f orini@physik . hu-berlin . de, f orini@pg . infn . it 



ABSTRACT: In the M = 4 super Yang-Mills theory, we consider the higher order anoma- 
lous dimensions 7l(<?) of purely gluonic operators TrT L where T is a component of the 
self-dual field strength. We propose compact closed expressions depending parametri- 
cally on L that reproduce the prediction of Bethe Ansatz equations up to five loop order, 
including transcendental dressing corrections. The size dependence follows a simple pat- 
tern as the perturbative order is increased and suggests hidden relations for these special 
operators. 
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1. Introduction 

Integrable structures emerge as a deep property of four dimensional Yang-Mills theories 
in the 't Hooft planar limit. In the simplest context, integrability underlies and governs 
the scale evolution of renormalized composite operators belonging to specific subsectors 
of the theory iJTJ] _ 

Historically, this intriguing phenomenon was discovered in the study of planar QCD, 
definitely a non-trivial quantum theory At one-loop, suitable maximal helicity Wilson 
operators admit a peculiar renormalization mixing matrix, the dilatation operator. It can 
be identified with the Hamiltonian of integrable XXX spin chains with s((2, M) symmetry. 
This is a light-cone subalgebra of the full four dimensional conformal algebra so(4, 2). 

^From a modern perspective, conformal symmetry, unbroken in QCD at one-loop, 
does not appear to be a necessary condition for integrability, as discussed in [^, ||, |5|, |6]]. 
Nevertheless, it plays an important role by imposing selection rules and multiplet struc- 
tures and is helpful to clarify the origin and details of integrability. The same reasoning 
applies to supersymmetric extensions of QCD with J\f = 1,2,4 supercharges. In particu- 
lar, multiplets of composite operators are greatly simplified in the maximal superconfor- 
mal M = 4 theory 0. Also, intermediate level integrability is achieved in various orbifold 
reductions of M = 4 SYM reducing the number of supercharges [Bj. 
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As is well known, in the maximal J\f = 4 case another conceptual tool is available 
to deepen the investigation, namely Maldacena AdS/CFT duality [^] . It relates M = 
4 SYM and AdS^ x S 5 superstring which is classically integrable [fC|]. Currently, a lot 
is known about the duality between the integrability properties of the two sides of the 
correspondence with a continuous very stimulating back and forth feeding. In particular, 
AdS / CFT duality has been a crucial ingredient to arrive at the higher loop proposal for 
the S-matrix of N = 4 SYM [glUl^, ^ |l^ [l^ g8[, |l^, glj, g^, g^, g4p. 

Forgetting for a while the string side, we can ask what kind of understanding can 
be gained from integrability in the gauge theory. This is a natural issue if we are ulti- 
mately interested in low-energy physical applications to hadronic phenomenology. As a 
first step, we can honestly postpone the important problem of taming conformal and su- 
persymmetry breaking and the precise link with QCD [P5fl. Working within the M = 4 
SYM theory, we can examine the outcomes and limitations of integrability as far as it is 
currently understood. 

,;From this point of view, integrability can be regarded as a tool for multi-loop com- 
putations, although this attitude could be admittedly narrow-minded. The typical object 
that is computed are higher order corrections to the anomalous dimension of specific com- 
posite operators. In all cases, these operators are single traces of the general form 

O = Tr ^Jl D nt Xij + permutations, (1.1) 

where X % are elementary fields in certain subsectors of the full M = 4 SYM and covariant 
derivatives generically appear to close the renormalization mixing. 

In the most favorable cases, we are able to write down Bethe Ansatz equations pro- 
viding the anomalous dimension of O as a perturbative series in the 't Hooft coupling g 

7o(5) = E c -(°)5 2n (1-2) 

n>0 

Such formidable results face a first and major limitation, namely the well-known wrapping 
problem (see [26, 27] for recent developments). The coefficients c„ are reliable up to a max- 
imum order n max that typically depends linearly on L. This means that 70(5) is actually 
calculable up to, say, 0(g 2L ) terms - a stumbling wall to any extrapolation to the genuine 
strong coupling regime. A notable exception occurs in the L — > 00 thermodynamical limit. 



Then, wrapping is absent and resummations of Eq. (1.2) can be attempted to match string 
duality predictions [|24], [2S| ]. 

A more subtle limitation appears when we try to investigate the dependence on L and 
{n t } at fixed perturbative order. Apart from very special cases, the Bethe equations do not 
provide the expansion coefficients as functions of L and {n;}, but just provide sequences 
of numerical (sometimes rational) values for each given operator. This is an unwanted 
situation as can be appreciated in the sector of the so-called twist operators [g^]. These 
are operators with a certain phenomenological origin in the QCD case. The length L is 



1 For those aspects of the duality that most concern our analysis, we refer the reader to [[10| fll| 
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fixed and one would like to know the analytic dependence of 7(5, S) on the spin quantum 
number S = J2 Z For instance, this is a standard procedure to analyze BFKL physics of 



pomeron exchange p0| ] 



Very recently intense work on twist-2 and 3 operators has led to higher order conjec- 
tures for the functions c n (S) appearing in the expansion 7(5, S) = J2 n c n {S) g 2n [^6[ [HI]. 
Proofs are however missing, at least beyond one-loop. It seems that new tools are needed 
to derive them rigorously from the Bethe Ansatz equations. 

If we give up exact results and turn to approximation, systematic methods can be 
applied to extract the large L, n; corrections. Indeed, this is a thermodynamical limit of 
the underlying spin chain where both the length and the number of magnons grows to 
infinity. Various techniques are available and have been successfully applied to rank-1 
5u(2) and s((2) subsectors ^ 0, g |5), |5eJ. For higher rank sectors (see also 



for an analysis in the rank-2 su(3) sector) the techniques developed in Q39| ] could be useful. 



Indeed, in the recent p0| ] an integral equation describing finite size corrections to the full 
nested Bethe ansatz was derived. 

On the string side of the AdS / CFT correspondence, there is an analogously intense 
ongoing discussion on how finite size effects of the string world-sheet could modify the 
solvability of the string sigma-model in AdSs x S 5 by means of a Bethe ansatz [Ml [13|, 
|T4| , f|2|, |3p. The currently known Bethe equations for quantum strings in AdS5 XS5 are 
asymptotic and describe the string spectrum with an exponential accuracy as long as the 
string length is sufficiently large [Jm|]. The breakdown of the asymptotic approximation 



via exponential terms, firstly described from a field theory point of view in [|45|], has been 
determined for the spectrum of spinning strings in the su(2) and s[(2) sector Q and 
for the giant magnon p7| ] dispersion relation in [|48|,^9|]. In particular, the exponential term 



in the finite size correction to the giant magnon dispersion relation has been recently and 



nicely rederived in [27] via a generalization of known results in relativistic quantum field 
theory, and there is a general and deep interest in obtaining exact results which should 
be valid for any value of the string lenght, which is in turn proportional to the lenght 
(i?-charge) of the corresponding gauge operator. 

In this paper, we contribute to the above general discussion and consider, in the gauge 
field theory context, a special class of operators where finite size corrections can be given 
in closed form. In other words, we provide the coefficients c n in Eq. fll.2| ) as exact functions 
of the operator length L. This is a seemingly unique result which, although peculiar, is 
very interesting and puzzling and certainly deserves some attention. 

The considered operators have a complicated mixing pattern and reduce at one loop 
to the purely gluonic higher dimensional condensates of the form 

O l = TtF l , (1.3) 

where T is one component of the self-dual Yang-Mills field strength. The operators Ol 
are exact eigenstates of the one-loop dilatation operator and can be mapped to the fer- 



romagnetic states of an integrable spin 5=1 chain |50|, plQ . As such, their one-loop 



dependence on the length L is trivial and (including the classical dimension), it is known 
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that [H, H^H] 

7 o L = L(2 + 3g 2 + ---), (1.4) 

Beyond one-loop, the analysis of [|2|] provide efficient computational tools to derive the 
sequence {c n (L)} for any given L, although not parametrically. It turns out immediately 
that c n (L) is not linear in L as far as n > 2. So, starting at two-loops, non-trivial finite size 
corrections appear. 

In this paper, we analyze the sequences {c n (L)} at fixed n as L is varied up to large 
values. By a careful investigation of the (infinite precision) numerics, we conjecture and 
provide closed expressions for c n (L) valid up to 5 loops, including the transcendental 
terms coming from the 5-matrix dressing phase JP7|]. To give an example, the two-loop 
anomalous dimension takes the remarkably simple form 

•ft W = 2 L + 3 Li + (-£ + | ( _ 1)L2 U + 1 ) + • • • , (1.5) 

with exponentially suppressed corrections to the trivial linear scaling with L. We have 
been able to extend the above equation up to five loops. The detailed results will be 
illustrated in the main text. Here, we just anticipate the large L limit which reads 

^ = M9) + S 4 Kg L) e- L ^ 2 + 0(e- 2 L lo * 2 ), (1.6) 



where fo(g) has been computed in p8|]. The function h(z) is regular around z = and 
does not receive contributions from the dressing phase, at least up to five loops. The size 
corrections to the thermodynamical limit are thus characterized by a finite specific corre- 
lation length £ = 1/ log 2. In the final Section of the paper, we shall try to argue why this 
correlation length abruptly appears at two-loops breaking the trivial linear dependence 
on L. 



2. One-loop ferromagnetic multi-gluon operators in the chiral sector 

In this Section, we introduce the special class of M = 4 multi-gluon operators that we are 
going to analyze. For completeness, we also review their one-loop integrability properties 
and, in particular, the reduction of the mixing matrix to the Hamiltonian of an integrable 
XXXi chain. 

In the planar limit, the most general purely gluonic local gauge invariant operators 
are easily identified. They are single trace operators built with covariant derivatives of 
the field strength 

TxiD^Fw-D^FwJ. (2.1) 

The anomalous dimension matrix T and the would-be spin chain Hamiltonian H are re- 
lated by 

where Z is the renormalization matrix and g is the scaled 't Hooft coupling kept fixed in 
the planar limit N c — > oo. 
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The one-loop operator H takes the form of a nearest neighbor Hamiltonian conserving 
the length L in Eq. It can be written 



H = #„,„+!, (2.3) 



n=l 



where the link Hamiltonian i?„ n+1 acts on the fields at positions n and n + 1 and is 
independent on n. For the complete = 4 SYM theory, the elementary fields are included 
in the singleton multiplet V JlO| ] and the link Hamiltonian reads f?3| ] 



= 2 £ MJ) Pf=\ Hi) = £ i (2-4) 

j=0 rc=l 



J 

n " 



where P^ =4 is a projector onto the irreducible superconformal multiplets appearing in 
the decomposition of the two-site states V ® V. 

To restrict the analysis to purely gluonic operators it is convenient to adopt the confor- 



mal analysis exploited in the QCD reduction described in [|5(J, tO]] (see also [p4|]). We first 
split into irreducible components the field strength F^ transforming as (1, 0) ® (0, 1) un- 
der the so(3, 1) = su(2) ® su(2) Lorentz algebra. This is achieved by means of the 't Hooft 
symbols [ p5| ] projecting F^ v onto the self-dual (1, 0) and anti-self-dual (0, 1) components 

= r\% f A + V$J A , A =1,2, 3. (2.5) 

The purely chiral gluon operators are the subset of Eq. dO| ) built using only the self-dual 
part of F,j, v 

Tr | D ni f Al ■ ■ ■ D nL J Al I . (2.6) 

At one loop, they close under renormalization mixing. The relevant link Hamiltonian can 
be obtained by restriction of H^ =i . To this aim, it is convenient to organize the various 
covariant derivatives of f A in a M = conformal infinite dimensional multiplet 

W = {D"/}»>o. (2-7) 

Two-site states decompose in irreducible multiplets labeled by the conformal spin j ac- 
cording to 

oo 

V s ® V f = V/ 1 . (2.8) 

j=-2 

Also, the conformal splitting of the full J\f = 4 projector P^ =i turns out to involve the con- 
formal projector Pj_ 2 only. This leads to the following purely gluonic link Hamiltonian in 
the chiral sector 

oo 

H = 2 £ h(j + 2)pf f . (2.9) 

Finally, if we further restrict to operators without derivatives, one can prove that the only 
modules appearing in the r.h.s. of Eq. (|L8|) are those with j = —2, —1, plfl. To make 
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contact with the spin chain interpretation, we introduce spin 5=1 su(2) operators {S 1 } 
acting on the three components f A as 

(S l f) A = ie lAB f B . (2.10) 

Then, the modules P- with j = —2, —1, can be shown to be associated with the su(2) 
representations with S = 0, 1,2 respectively, appearing in the decomposition 1 ® 1 = 
2 ® 1 ® 0. The link Hamiltonian H n ^ n+ i in Eq. (^9|) can be written as a polynomial in 
S n ■ S„_i_i with the result 



H = 3L + \ E [ s « • s "+i - ( s « • s «+i) 2 ] • ( 2 - n ) 

n 

This is an anti-ferromagnetic integrable spin-chain that can be diagonalized by Bethe 



Ansatz [56|. The ground state is highly non- trivial, but the maximally excited states are a 



trivial ferromagnetic multiplet. A convenient representative is the operator 

O l = Tr{F L ), (2.12) 

where T = f + is the maximal eigenstate of S z . The anomalous dimension of this state 
(including the classical dimension) is simply 

lL (g) = 2L + 3Lg 2 + 0(g 2 ). (2.13) 

The linear dependence on L follows from the uniform structure of the ferromagnetic state. 

In this paper, we shall be working on the conformal/ field-theory side of the AdS / CFT 
correspondence. However, it must be mentioned that the natural candidate for a semi- 
classical string state dual to Ol has been proposed in [£7L[53]. It describes a rigid circular 
string rotating simultaneously in two orthogonal spatial planes of AdSc, with equal spins 
S = L. At large S, the weak-coupling extrapolation for the energy is given by 



E = p(A) S + q(X) + ..., 
Pi 

Vx 



Px>i = Po + -7= + • • • , 3a»i = v^go +?! + ... • (2.14) 



The results for the three-level and 1-loop coefficients of the solution have been calculated 
in p8[ ] within a stability region 0.4 < -j= > 1.17, corresponding to a fixed value (m = 1) of 
the winding number. It should be noticed that, since in the semiclassical approximation A 
is large on the string side, the interval of stability for the solution does include large values 
of S, allowing the comparison to large S, large A asymptotics of the exact anomalous 



dimension. The linear dependence on S exhibited by the solution (|2.14| ) supports the 
identification of the gauge theory operator Ol with this particular rigid spinning string 
solution. 

3. Higher loop extension of the scaling field O l 

At more than one-loop, the operator Ol ceases to be an eigenstate of the dilatation oper- 
ator because the purely gluonic chiral sector does not close under mixing anymore. The 
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higher order scaling operator receives corrections and contributions from the other sectors 
of the full p5u(2, 2|4) theory and is uniquely defined by the boundary condition of being 
Ol at one-loop. In the following, we shall not be pedantic about this distinction and keep 
naming Ol the multi-loop extension of TrT L . 

A great deal of information about Ol has been obtained in [[Z|] in the framework 
of the long-range Bethe Ansatz equations. In this Section, we quickly summarize these 
results with some additional investigation of the finite but large L Bethe roots. This will 
fix the setup for the computation of 7l(<7)- 

3.1 Dynkin diagrams and Bethe roots 

As is well known, several choices are available for the Dynkin diagram of a Lie superal- 
gebra. In the case of psu(2, 2|4), the one loop analysis of the operators Ol is almost trivial 
with the Kac distinguished form. Indeed, Ol is the vacuum state and no calculation is 
needed. 

On the other hand, the all-loop Bethe equations are known for a limited set of (differ- 
ent) choices of the Dynkin diagram [[l6|]. In particular, we shall work with the following 
one 

+1 

—(g) O <gn^>Hg> O <g)— O.I) 

With respect to this Dynkin diagram, the vacuum is the BPS state Tr Z L and Ol is a highly 
excited state with many excitations, whose momenta have to be diagonalized by solving 
the Bethe Ansatz equations in order to reproduce the correct energy. The excitation pattern 
of Bethe roots for Ol is 

{K U K 2 , K 3 , K i} K 5 , K 6 , K 7 ) = (0, 0, 2L - 3, 2L - 2, L - 1, L - 2, L - 3) (3.2) 

where Ki is the excitation number of the i-th node of the Dynkin diagram 

+1 

—(g) O <8^>ng> O ®— 0.3) 

2L-3 2L-2 L-l L-2 L-3 

All but the first two nodes are highly excited. 

All the one-loop Bethe equations can be exhibited as roots of explicit polynomial by 
means of the dualization procedure illustrated in |5^] to which we defer the reader for 
more details. It is instructive to describe the procedure in graphical terms. Dualizing first 
at nodes 3 and 7, we obtain 

+1 

— ® <g^Hg> <gHHg> <gHHg) (3.4) 

2L-2 L-l L-2 
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Dualizing at nodes 4 and 6, we obtain 



+2 



L — 1 



O 



(3.5) 



The Bethe equations are thus reduced to the simple equation 

L 

I =1, k = l,...,L-l, 



U 5 ,k + i 



U5,k 



which is solved by 



U 5> k = cot 



7T k 



(3.6) 



(3.7) 



The dualization process can be inverted step by step providing exact polynomials whose 
roots are the Bethe roots at any finite L. In particular, one finds for the roots at node 4 and 



6 the explicit result p8| ] 



Qi{u) 



u + 



u + 



3i 



u 



+ u 



3i 



u 



u + 



3i\ L ( 3i xL 



u 

2 



(3.8) 



Of course, from the knowledge of Qi(u) one can prove again the one-loop result Eq. fl2.13| ). 

Going over to higher orders, we have to work with the long range Bethe equations 
which are a deformation of the one-loop ones. They involve the standard quantities 



x(u) = — 



1 + 



2g 2 



X ( u ± - 



(3.9) 



and read 



C 4,fc, 



2L-2 

n 



^3, fc - 



2L-2 

n x 

j^k 
2L-3 



: t,k X i,j 1 9 2 /^ x t,k X i,j 

x t,j 1 ~~ 5 2 /2 x ^k x t,j 



: 4,fc X 3,j 



L-l 

n 



4,fe x 5,j -p-|- 1 g /2x A k X7j 
g 2 /2x+x 7 /' 



3J j=l 1 



j=l :E 4,fc X 3,j J=1 ^4 ( fc 2l5,i 

£—2 „, „, i i 2L—2 _ _+ 

j-j" u 5,k - Uqj + | 3=5, fc - a 4| j 

j =1 U 5ik - U 6 ,j - \ j =1 %5,k ~ 

j-\ ^6,fc - U6,j - i rr u 6,k ~ 1*5, j + f pi- 3 ^6,fc ~ U 7 j + | 

3 - j u 6,fc - U 6f j + 1 ^ U 6 , fc - U B j - I j} x u 6tk - U 7 j - \ 

L-2 

n 



(3.10) 
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where a 2 (uk,uj) is the dressing phase to be discussed later. 

It is possible to perform a partial dualization of these equations and obtain reduced 
long-range equations involving roots at nodes 4, 5, 6 only. These are 

C, \ L 2.Zj 2 ~~\~ 1 — 2Zj 4 ~\~ > 

S 4,fc 1 _ I"! 3: 4,fc ~ a: 4,j 2x Zk x 4,i J2, ni \ T-f a! 4,fe ~ X 5J 

^.fc/ j=l ^.fc a3 4 ) j 1 ~ _+ i=l JC 4,fe X 5,j 
jjtk 4,fc 4,i 

i = n ! 5 - fc = U6j + _ \ n 2 ! M ~ a y , o.id 

j=l u 5,k u 6,j 2 j=l X5 < k X i,j 

1 _ T-r U%,k - U 6 j + I -|-r ^6,fc - Ub,j ~ 2 
/J[ ^6,fc - u 6j - i ^ u 6ifc - U S J + § ' 

Here, S5 are the 2L — 4 roots dual to U5 . At one loop, they are the roots of the polynomial 

Q B (u) = 3u 2i + (-< + u) L {-2i + u) L + (2i + u) L ({i + u) L + (-2i + u) L ) 

((-i + u) L + (t + u) L + 2(-2i + u) L + 2(2t + u) L ) . (3.12) 



The roots at nodes 4 and 6 are still given at one-loop by the previous polynomials. At 
generic g, the anomalous dimension is obtained from the roots ^4^(3) alone and reads 

nig) = 2 l + g 2 f: ( - —^) ■ (3.13) 

Finally, let us consider the dressing phase. It enters the calculation starting from four 
loops. Its general form is discussed in [^4j]. The terms relevant for a computation up to 
five loops are simply 

a 2 {u,u') = e^ { - u ' u '\ (3.14) 

where 

ti(u, u') = (C 3 3 6 - 5 Cb 3 8 ) (? 2 (tt) <7sK) - 92K) g 3 (tt)) + • • • , (3.15) 



v x+(u) a;-^)/' 2 U + («) 2 x-{u) 2 J ' 

The coefficients are transcendental sums £ n = £~ n - 

3.2 The one-loop Bethe roots: some numerics at large but finite L 

The one-loop Bethe roots are the zeroes of the polynomials Qi$fi(u). It is instructive to 
study them at large L comparing with the results of [|2^] obtained in the L — » 00 limit. First 
the (dual) roots u s . They are complex. We show them at L = 100, 200, 350, 500 in Fig. (@). 
As predicted, most of them are distributed along two segments with Im« 5 = ±|. Apart 



-9- 



from these roots, other ones are scattered in the complex plane according to a nice regular 
pattern. To understand these roots, we look for a Bethe root admitting the expansion 

u = L (x + ^ + ^ + ---y (3.16) 

A roots with leading behavior u ~ L is called extremal in [^8|| . Replacing this expansion in 
Qe (u) , we obtain a well defined large L expansion for the ratio 

R( X0}Xl ,...;L) = ^^-. (3.17) 

The leading term is 

R = 64 cos 2 — sin 4 — + ©(L" 1 / 2 ), (3.18) 
2xo 2 so 

leading to xq = -p- for any integer n. Considering separately the cases n even/ odd and 
expanding at higher order in L -1 / 2 one finds the solutions (in the first quadrant) 



a± ' k ~ (2fc + l)7r 



L±J-\+o(L- 1 / 2 ) i A; = 0, 1,2 (3.19) 




P±,k = ^ (L±±VZ^3±iVlS±-^) +o(L-V 2 ), k = 1,2,.... (3.20) 



The other roots are related by reflection with respect to the coordinate axis. For large L, 
the a-roots appear in real close pairs. These pairs are closer to the origin as k is increased. 
In general, for a given L, only a finite number of such pairs is well approximated by the 
above formula. The /3-roots have an imaginary part and also appear in close pairs. In 
Fig. (||) we draw crosses at the first a and /3 pairs. 

The roots u^ in and Ue t n are real. Their density is defined in the L — > oo continuum 
limit as p(u) = dn/du and the analytical prediction is 

Pi{u) = h (^TI + ^Tf) ' P6iu) = ^^TJ- (3 - 21) 



In the discrete case at finite L, we can plot the points 

'u n + u n+ i 4L 



u n+l ~ u r 



(3.22) 



The result is shown in Figs. for the Bethe roots at L = 200. The agreement is quite 
good in the case of u^. For u%, we observe a deviation in the tails of the distribution at 
large |ua|. It can be understood as in the above discussion of extremal U5 roots. 

4. Perturbative expansion of the long-range Bethe equations 

Starting from the exact (i.e. known with arbitrarily high precision) one-loop Bethe roots 
we can make a perturbative expansion in even powers of g 



u a , k = J29 2n u%l a = 4, 5, 6, k = 1, . . . , K a , (4.1) 

71 = 
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where we relabel U5 = U5, K$ = K$. The explicit five loop expansion of the anomalous 
dimension can be compared with the results of [^] up to L = 8. We have extended the 
calculation up to L = 60. The results at five loops for L < 20 are shown in Appendix (A). 

The expansion is a rational combination of 1, £3 and £5. As we mentioned, the zero 
and one loop results are proportional to L. In general, it is convenient to redefine 

lL {g) = L(2 + Zg 2 + Y J Cn{L)g 2n \. (4.2) 

We now show that it is possible to provide simple closed expressions for the non-trivial 
functions c„(L). As a constraint, we must meet the exact expansion in the L — > 00 limit 
obtained in p8| ] and reading at five loops 

c 2 (oo) = -—, 

c 3 (oo) = (4.3) 
, . 59487 27 , 

Ca(OO) = C3, 

4v ; 512 4 sd, 

, . 632661 1665 , 135 , 

CrIOO) = 1 C3 H C5- 

bv J 1Q24 32 sd 4 s& 

As a general remark, it is instructive to plot the numerical values of c n (L) at the first values 
of L. Indeed, it is immediately clear that factors (— 1) L can appear in the closed formula 
for c n (L). Therefore, we shall analyze the odd and even L cases separately. 

4.1 Two loops 

For odd L = 5,7,9, ... , we subtract the asymptotic value 02(00) and rescale to find 

s (L) _ c (oo)) = _ 1 _1_ _ J 1 1 1 1 l_ 

3 ( i c 2 (, j c 2 (oo)) g> 2i , g5 , 34i> i3g5 , 54gi> 21g45 , g7381 . 349525'- " 

(4.4) 

A careful inspection reveals that the denominators are simply related to powers of 2 minus 
one. The precise formula is easily found and reads 

c 2 (L) = -—-- —i , L odd. (4.5) 

8 8 2 1 -- 1 - 1' 

We checked it for all the L that we have explored. Remarkably, it works also for the even 
L case if the sign of the term ~ 2 L is changed. The final formula is thus 

^) = -T + i (-i^2^ + i - (46) 

This simple result is rather remarkable. It holds at finite L and predict exponentially 
suppressed deviations from the trivial linear scaling of the anomalous dimension 7 ~ L, 
valid up to the one-loop level. Is it possible to obtain a similar result for the next three 
loop contribution ? 
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4.2 Three loops 

Following the strategy adopted in the two-loop case, we start again from odd L = 5, 7, 9, . . . 
and evaluate 

, \ HI 425 3628101 _ 

c 3 (L) - c 3 (oo) = H , . , (4.7) 

v ; K ' 1600' 3136' 39304000' 

9904623 7874523 63804855621 



(2 L_1 - l) 3 (c 3 (L) - c 3 (oo)) = 185, -26775, -1209367, -36316951, -921319191, (4.8) 



230701504' 463736000' 10423090379584' 
This sequence appears to be definitely non trivial and much more complicated than the 
two-loop case. In particular, the signs are not definite and the denominators do not have 
simple factorization properties. However, the sequence enjoys a remarkable property. If 
we multiply it by (2 i_1 — l) 3 and apply a constant scaling, we find 

2_ 6 
3 4 

-21268285207, -461958727447, -9613145655063, -193758643734295, .... 

Indeed, the sequence is integer for all considered L. As a second feature, one can plot the 
following function of L 

(2 i - 1 -l)(c 3 (L)-c 3 (oo)), (4.9) 

and it turns out to be curve quite close to a quadratic parabola. From these two features, 
it is natural to look for a closed formula of the form 

2 2 

ca(L) - c 3 (oo) = I E 2PL E C ™ L<? - ( 4 - 10 ) 

Indeed, it turns out that all the three loop results at odd L are reproduced by 

, r . 393 -9 • 2 2i (9L 2 - 33 L - 104) - 18 • 2 L (9L 2 + 15L + 202) + 3528 r 

c 3 (L) = 1 \ , L odd. 

V ' 16 64(-2 + 2 L ) 3 

(4.11) 

Looking back at Eq. fl4.6| ), there is a striking similarity suggesting an all order structure. In 
particular, the same formula works for even L, if we apply the modification rules 

2 2pL 2 2pL , 2 (2p+1)L -> -(-l) L 2 (2p+1)i . (4-12) 

The general formula is then 

393 9 • 2 2L (9 L 2 - 33 L - 104) - 18 • (-1) L 2 L (9 L 2 + 15 L + 202) - 3528 
c 3 ( J-^g"+ 64-8 [(-l) L 2 i " 1 + l)] 3 ' 

(4.13) 

4.3 Four loops 

At four loops, we attempt to repeat the game. The only new feature is the transcendental 
contribution from the dressing phase. This is a piece of c$(L) proportional to £ 3 . From the 
numerics, it is independent on L and reads 

c 4 (L) = 4 0) (L) + ci 3) (L)C 3l 4 3) (L) = - (4.14) 
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The c^{L) is a rational contribution with properties quite analogous to those of c${L). In 
particular, for odd L 

(0 ^(2 i - 1 -l) 5 (cl 0) (L)-4 0) (oo))eN I (4.15) 
(ii) [2 L ~ 1 - 1) (c 4 0) (L) - c 4 0) (oo)) ~ L 4 , for L — > oo. (4.16) 
Again, it is natural to postulate from (i) and (ii) the closed formula 

4 0) (l) - 4 0) (oo) = ? - 3 * I E 2?L E ( 4 - 17 ) 

Replacing the explicit anomalous dimensions in this formula we find that indeed it is 
satisfied by all considered (odd) L with coefficients d p<q giving 



(0) 59487 2 4L g 4 , 4 + 2 3L Q 4 ,3 + 2^g 4 , 2 + 2^4,1 - 1335168 
" 4 (L) " "~512~ + 2 i5 (2 l-i _ 1} 5 ' L ° dd ' (418) 

where the Q polynomials are 

g 4 ,i = -72 (27 L 4 + 90 L 3 - 1485 L 2 - 2004 L - 38456) , 
g 4 , 2 = -108 (99 L 4 - 18 L 3 + 513 L 2 + 2958 L + 19924) , (4.19) 
g 4i3 = -18 (297 L 4 - 1080 L 3 + 1647 L 2 - 12060 L - 41264 
g 4 , 4 = -9 (27 L 4 - 252 L 3 - 1053 L 2 + 5190 L + 10676 



The case L even is obtained changing the sign of 2 L in the powers (2 L ) P and correcting 
with a shift in the boundary case L = 4. The final result is 

(0 ), M 5 59487 , 

C 4 ' L = dr, 4 h 

4 v ; 64 ' 512 

2 4L g 4 , 4 - (-i) L 2 3L g 4|3 + 2 2L g 4 , 2 - (-i) L 2 L g 4|1 - 1335168 

2 15[(_ 1 )L 2 L-l + 1 ]5 

The above shift, as well as the corrections appearing in the five loop formula ([L26) be- 
low, are possibly related to short wrapping effects - the lack of the asymptotic conditions 
prevents in the boundary cases the validity of the Bethe equations. 

4.4 Five loops 

At five loops, we have a more complicated dressing contribution with two different tran- 
scendentality terms 

c 5 (L) = 4 0) (L) + 4 3) (L) C 3 + 4 5) (L) C 5 , 4°' 3,5) (L) e Q. (4.21) 
The maximum transcendentality is independent on L 

cf'W = ^. (4.22) 
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Repeating the above heuristic analysis for the other terms we find for the transcendental- 
ity 3 term 

4 3) W = ^ - I h * + 81 -(-1)^-4) - 648 
The purely rational term has the representation 

4 0) (l) - 4 0) (oo) = 7 - r * e 2?i E e ™ Lq > ( 4 - 24 ) 

with the explicit final result, holding for even or odd L 

(0) _ 632661 14987 333 G(L) 

Cb " ~T024~ + 12288 ^ L ' 4 " 4096 8l >* + 2 ™ [(-l)i2^ + lf ( ) 

G{L) = J2(-lY {L+1) 2 pL Q5,r, (4-26) 

p=0,6 

where 

Q 50 = -2 11 • 3 2 (432 L + 56639), (4.27) 
Q 5i i = 2 5 • 3 4 (9 Z, 6 + 45 L 5 - 1581 L 4 - 4113 L 3 + 39492 1? + 53316 L + 1253696) , 

Q 5)2 = 2 4 • 3 3 (l539 L 6 + 1755 L 5 - 37503 L 4 + 41409 L 3 - 370980 1? - 961116 L - 9751792) , 

Q 5)3 = 2 4 • 3 2 (12231 L 6 - 17172 L 5 + 68067 L 4 + 158976 L 3 + 358722 Z, 2 + 3589416 L + 20219128) , 

Q 5 , 4 = 2 3 • 3 5 (453 L 6 - 1956 L 5 + 3769 L 4 - 15096 L 3 + 13278 1? - 163616 L - 582008) , 

Q 5)5 = 2 • 3 3 (l539 L 6 - 13095 L 5 - 10611 L 4 + 82683 L 3 - 290952 Z, 2 + 1783716 L + 4340656) , 

Q 5 , 6 = 9 (81 L 6 - 1377 L 5 - 6129 L 4 + 103653 L 3 + 195912 L 2 - 1277388 L - 2247232^ 



The extension to higher loops seems to be a computational issue. One has to generate 
a large enough number of terms in c n (L) and must check that an Ansatz similar to the 
previous ones matches it. 

5. Large L expansion of 7l(s0 

The five-loop results described in the previous sections are valid at finite L. Nevertheless, 
it is interesting to look at the dominant terms at large L. As remarked in the Introduc- 
tion, the resulting expression can admit a thermodynamical interpretation. Collecting the 
formulae for c n and expanding at large L, we find 

^ = /o(s) + fi(9, L) e~ L lo * 2 + f 2 (g, L) e~ 2L ^ 2 + f 3 (g, L) e~ 3i lo s 2 + • • . . (5.1) 



The leading term agrees by construction with the result of [28] 



* , x „ 2 51 4 393 6 / 27Cs 59487\ 8 /1665Cs 135& 632661 \ 10 

fjq) = 2+3 q 2 o 4 H o 6 + — o 8 + — H — H q w +. . . . 

joyyj y s y ^ iq y ^\ 4 512 ^ * V 32 4 1024 J y 

(5.2) 
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The first exponentially suppressed term has a prefactor 

. . r , 9 4 / 81Z, 2 297L 117\ 6 
M9,L) = - l9 *+[-— + — + — j 5 6 + 

243L 4 567L 3 9477L 2 23355L 24021 \ 



+ 1 1024 + 256 + 1024 512 256 



" 9 + 



( 729L 6 12393L 5 55161L 4 932877L 3 220401L 2 
+ 1 ~ 32768 + 32768 + 32768 32768 409(5 + 



+ 



/2874123 _ 8IC3 X 81& 316017A 

V 8192 32 ) 8 512 / y v ' 



At order 0(g 2n ), the leading power of the length is L 2 "~ 4 and comes always in transcen- 
dentality terms unrelated to dressing. The large L limit of fi(g, L) can be compactly 
written as 

/i(5,i) = -^ 4 (l + ^ 2 + y + ^ + •••), * = \l9, M 

and in particular, given the absence of transcendental contributions, do not depend on 
the dressing phase. It seems reasonable that this structural properties could persist at all 
orders. 



6. Discussion and Conclusions 

In this paper, we have considered the chiral operator Tr T L in M = 4 SYM. At one-loop, 
it scales with a definite anomalous dimension j L proportional to L. At two-loops and 
beyond, it mixes with the other psu(2, 2|4) fields. The length L is no more a conserved 
quantity and 7l /L is not constant. In principle, this ratio is not expected to be expressed 
by a simple expression at finite L. One would just resort to compute systematically its 
corrections at large L. 

Nevertheless, the main result of this paper shows that some unexpected structure ex- 
ists at finite L. We have been able to provide a closed form for Jl/L up to five-loops. 
Radiative corrections follow a simple pattern order by order in perturbation theory, in- 
cluding transcendental dressing effects. They are sensitive to the parity of L and are ex- 
ponentially suppressed as L — > oo. 

A closed formula for the multi-loop size dependence is a remarkable fact that has 
no counterpart in existing calculations for other operators in the various subsectors of 
M = 4 SYM. It can be due to the simplicity of the considered operator or could hint to 
some hidden relation obeyed by the anomalous dimensions as a function of L. The closed 
formulae are a mere conjecture, although with a strong empirical basis. It is clear that a 
(dis)proof would be certainly enlightening. 

In the large volume regime our results read 

^ = fo(9) + g*h{g L) e- L '°e 2 + • • • . (6.1) 
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Eq. claims that starting at two-loops, exponentially suppressed corrections appear 
with a g independent correlation length £ = 1/ log 2 and the combination g L as a nat- 
ural scaling variable for the prefactor. It would be interesting to understand such fea- 
tures from the point of view of the spin-chain interpretation of the dilatation operator H. 
We emphasize that the C(2~ L ) corrections have nothing to do with much smaller C(A L ) 
wrapping effects. A natural explanation for the exponential corrections could take into 
account length-changing processes as suggested in p8|]. An explicit two-loop calculation 
of H would be important to clarify these issues. 

We conclude with a remark concerning the dressing phase i?. Currently this is a 
well understood ingredient appearing in the 5-matrix. However, it would be very nice to 
classify the special kind of interactions that are associated with it in the dilatation operator. 



A relevant step in this direction has been recently described in [60] where it is linked to 
so-called maximal reshuffling interactions. In our investigation, the special feature of 
dressing effects is that they are subleading at large L and up to five-loops. Transcendental 
contributions drop out from the function fi(L,g) being characterized by subdominant 
powers of the length L. 
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A. Five loop anomalous dimensions for L < 20 

/ x o ,„ 2 or 4 1515 6 ( 513937 n „ \ 8 
74(5) = 8 + 12 5 2 - 25 5 4 + — 5 6 + 27 Cs J $ 8 

/ 22129823 1651 inr ^A in 

' Ki)> * 4 320 V 12800 4 V 

/ 1539949881 8307 675 \ 10 

+ l 512000 + ^2" 6 + ^ Cb J 9 

2 837 4 6278355 fi / 29266837713 81 , 

76 (0) = 12 + 18 g 2 g 4 H g 6 + & 

myj y 22 y 42592 y V 41229056 2 w 

/ 76857234976107 605205 405 \ 10 

+ V 19954863104 + 1936 " + ~2~ 5 + " ' ' 

9 179 4 76603 fi / 181131695 189 \ 8 

77(g) = 14 + 21 g 2 g 4 + g 6 + 6 5 

uwi v 4 448 V 225792 4 / 

Z8959397257 11641 945 \ 10 

+ V 2107392 + ^T C3 + — (5 ) 9 + -"' 

2 2190 4 125533809 fi / 8840715968859 

7 8 (g) = 16 + 24 g 2 g 4 H g 6 + 54 6 

,8vyv y 43 y 636056 y \ 9408540352 v 

/ 346753221469919673 3080871 2 \ 10 
+ V 69585564443392 + 7396 ^ 3 + ^ J 9 + " ' ' 

„ 2 19521 4 8655987591 fi / 2364798793587021 243 , « 8 

79(g) = 18 + 27 g 2 g 4 H g 6 + 6 g 

' vyy * 340 y 39304000 y V 2271771200000 4 3 ' b 

( 73033337654861466627 108214623 1215 \ 10 

+ V 13130837536000000 + 231200 ^ 3 + 4 ^ 5 J 5 + " ' > 

2 7265 4 486455845 fi / 179336215108445 135 , 

7io(5) = 20 + 30 g 2 g 4 H g 6 + (3 

' y * J y 114 * 1975392 y \ 154033166592 2 s 

/ 4103422374381475165 27055595 675 \ 10 

+ V 667271677676544 + 51984 " + ~2~ ^ 5 ) 9 + " ' ' 

2 8697 4 5656701069 fi / 206094402320199 297 \ 8 

711(g) = 22 + 33 g 2 g 4 + g 6 + £3) 9 

' Kyj y 124 y 20972864 y V 161239378432 4 s / * 

/512022650767888797357 17597781 1485 \ 10 

+ V 74996304653805568 + 30752 ^ 3 + 4 ^ 5 J 5 + " ' > 

2 52245 4 1504241745363 fi / 26503161491873431953 

712(5) = 24 + 36 g 2 g 4 H g 6 + 81 6 

' 683 y 5097791792 y V 19024510362066304 St 

/ 262052003573439955673753835 2330333685 405 /^ 10 

+ V 35498899257159792346624 + 3731912 ^ 3 + ^ 5 ) 9 + " ' ' 

2 H603 4 11382640977 fi / 61846229508401447 351 , 

713(g) = 26 + 39 g 2 g 4 H g 6 + 6 

' Kyj y 140 y 35672000 y V 40901515200000 4 s 

M197004793411747026501 26513707 1755 \ 10 
+ V 521085303648000000 + 39200 ^ 3 + 4 5 + " ' ' 

2 487473 4 224231872961943 fi 

714(g) = 28 + 42 g 2 g 4 + g 6 

' K * J 5462 y 651801084512 y 

/ 63195216734569435123965 189 
+ V 38890942307856038656 2~ ^ 3 J 9 

( 10018082258604066826691362529007 86929777809 945 \ 10 
+ 1 1160250749448653889305611264 + 119333776 ^ 3 + ~Y ^ 5 J 5 + " ' 



8 
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2 2088855 4 3839300288893665 fl 

715(9) = 30 + 45 g 2 g 4 H g 6 

/ibw/; -r y 21g44 y -r 1042 3090379584 y 

/ 4337543080186955113152555 405 \ 8 

+ V 2486742653840334490112 4~ " ) 9 

/2751363092765611488400885340297145 744805296585 2025 \ 10 

+ V 296643740062996423913557649408 + 954320672 " + 4 ^ 5 J 9 + " ' 
„ 371380 . 75888854970083 „ 



2 371380 4 75888854970083 fi 

7i6 (5) = 32 + 48 g 2 g 4 H g 6 

' * 3641 y 193073214884 * 

/ 342362914175507036271937 , \ , 

+ 108 6 g' 

V 184287501648332409888 s J * 

( 10731015336070384155651542820939 22072903471 \ 10 

+ V 1085812212950776269479306368 + 26513762 ^ 3 + ^ 5 ) 9 + " ' ' 

2 557049 4 963880152452127 fi 

717(0) = 34 + 51 g 2 g 4 + g 6 

' w * 5140 y 2308544648000 y 

/ 1024443318894708985760901 459 \ 8 
+ V 518422022549556800000 4~ " J 5 

/ 122265300022822344725447154392091 46737698103 2295 \ 10 
+ V 11642027096907730208288000000 + 52839200 ^ 3 + 4 5 + " 

2 10027071 4 1180027266090140025 fi 

7l8 (g) = 36 + 54 g 2 g 4 H g 6 

' K * J * 87382 * 2668860863627872 y 

/ 85208553177354546108032208051 243 \ 8 
+ V 40756782343071289195379456 2 ~ " J 5 

/3461778112225195512629525977023541804641 28604992980447 
+ V 311203054756192480704870130697145344 + 30542455696 ^ 3 

\ 1 n 



1215 f x 10 



2 / 

„ 2 742739 4 681479838093317 fi 

719(g) = 38 + 57 g 2 g 4 H g 6 

' * 6132 y 1460288902464 * 

/ 3455324577868750102309487 513 \ 8 
+ V 1564904852245239998976 4~ " J 9 

( 1093345538309226287514039899518949 74344962037 2565 \ 10 
+ V 93167530542473483540155557888 + 75202848 " + 4 5 + " ' > 

2 22282275 4 41954422215573649845 fi 

720(5) = 40 + 60 g 2 g 4 H g b 

' w 174763 * 85402081606607152 y 

/ 48481600634126052932187928286535 , \ 8 

+ 135 6 3 

V 20866875547860781826069765504 s J y 

/ 31509537567806883433426191930418295740388445 254263222600515 
+ V 2549273313592296157452673880096407176704 + 244336849352 ^ 3 
+675 Cb) 3 10 + -- - , 

(A.l) 
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Figure 1: Dual Bethe roots u 5 computed at one-loop with L = 100, 200, 350, 500. Crosses on the x 
axis are pairs of a± extremal roots. Crosses with non zero imaginary part are j3± roots. 
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Figure 2: Density of Bethe roots u 4 from the analytical prediction p 4 and from the numerical roots 
at L = 200. 




Figure 3: Density of Bethe roots u 6 from the analytical prediction p 6 and from the numerical roots 
at L = 200. 



